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STATISTICS AND MULTIWAVELENGTH SYNTHESIS MODELS: 
TOWARDS A NEW GENERATION OF SYNTHESIS MODELS 

M. Cervino 1 ' 2 
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RESUMEN 

En esta contribution presento mi trabajo en el desarrollo de una nueva generation de modclos de smtesis 
evolutiva que calculan tanto el espectro multirango (desde rayos-7 a ondas de radio) como la dispersion asociada 
en cumulos estelares jovenes. Tambien muestro algunos efectos estadfsticos que pueden aparecer en en analisis de 
conjuntos de datos, como distribuciones bimodales o multimodales y sesgos cuando se comparan los resultados 
de los modelos de smtesis con datos reales. Esta nueva generation de modelos de smtesis puede ser una 
herramienta muy util para el analisis de los datos que se esperan obtener con el GTC. 

ABSTRACT 

In this contribution I present my current work in a new generation of evolutionary synthesis models that com- 
pute the multiwavelength energy distribution (from gamma-rays to radio) as well as the associated dispersion 
for young stellar systems. I will also show some statistical effects that may appear in the analysis of surveys, 
like bimodal or multi-modal distributions and bias when color indices computed by the codes are compared 
with observations. Such new generation of synthesis models may be useful for the analysis of the data expected 
from GTC. 
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1. INTRODUCTION 

Taking into account the discreteness of the stellar 
population, the predictions of any model that relies 
on an Initial Mass Function (IMF) are only exact 
under the assumption of an infinite number of stars. 
Otherwise, they only give a mean value of a probabil- 
ity distribution. The relevance of such fluctuations 
in the results of synthesis models is obvious in the 
case of massive stars and young clusters, but they 
also affect the models of older clusters dominated by 
the emission of low-mass stars since small variations 
in the initial mass/number of stars in a given mass 
range, can produce different numbers of, e.g., AGB 
stars at a given age, which in turn produce large 
variations in the resulting observable. 

So, for the comparison of models with observa- 
tional data it is necessary to obtain not only the 
mean value of the observable, but also, at least, the 
corresponding dispersion of the computed observable 
due to the discreteness of the stellar population, and, 
ideally, the underling probability density distribu- 
tions. The dispersion due to a finite number of stars 
in real systems can be evaluated theoretically, as it 
has been shown in Cervino et al. (2002). 
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2. SOME EXAMPLES: MULTIWAVELENGTH 
EMISSION AND COLOR-COLOR DIAGRAMS 

In Fig. [l] I show the 90% confidence level for the 
multiwavelength spectrum for a 5.5 Myr old burst 
obtained analytically from our code. Examples of 
statistical effects of a large number of Monte Carlo 
simulations and comparisons with the dispersion ob- 
tained theoretically for different observables can be 
found in Cervino et al. (2000, 2001, 2002). 

In Fig. H I show the B-V vs. V-K colors for a 10 
Myr old star forming region. Note that it is assumed 
that colors are independent on the number of stars 
in the cluster, so I have test this hypothesis with 10 3 
Monte Carlo simulation of clusters (triangles) with 
1000, 100, 10 and 1 stars in the mass range 2-120 
M Q following a Salpeter IMF slope. The theoreti- 
cal result is shown with the symbol. The Figure 
shows two important effects: (a) There are bimodal 
(multi-modal) distributions of the colors, (b) The 
mean value obtained by infinite populated IMF (an- 
alytical result) has a bias if it is used for the analysis 
clusters with a small number of stars. As reference, 
the mean mass of a cluster with 10 3 stars in the 
given mass limits is w 6 x 10 3 M Q , it means 2 x 10 4 
M Q if the value of lower mass limit of the IMF is 
0.08 M Q . Both, bimodality and bias, are a natural 
effect of statistics when only a few stars dominate 
the observable. Note that it is dependent on the ob- 
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Fig. 1. Analytical 90% Confidence Level for the multiwavelength spectrum for a 5.5 Myr old burst. 




Fig. 2. Monte Carlo simulations of the B-V vs. V- 
K colors of 10 Myr old stellar clusters with 1000, 100, 
10 and 1 stars in the mass range 2-120 Mq following a 
Salpeter IMF slope (triangles). The theoretical result is 
showed with the symbol. 

servable: Example 1: A cluster with 10 5 stars in the 
mass range 2-120 M© with a Salpeter IMF slope 
have about only 10 WR stars!. Example 2: The 
effect on the colors of globular clusters (dominated 
by a few AGB star) can be relevant (see Bruzual 
2001 and references therein for details). Addition- 
ally, these bimodality effects produce a bias when 
ratios and logarithm quantities computed by synthe- 
sis codes are compared with observed data (but such 
a bias is not present for quantities like luminosities 
even for "clusters" with only 1 star). This subject is 
addressed in Cerviho & Valls-Gabaud (submitted). 



3. DISCUSSION 
Statistical effects may be specially relevant (i) 
In the analysis of individual systems. In particular 
for the galaxy building blocks if they are small clus- 
ters, (ii) In the analysis of the integrated spectra 
obtained with high-resolution instruments when not 
all the region is covered by the slit. (Hi) In chemical 
evolutionary models (Cerviho & Molla, submitted), 
as far as they are dependent on the yields production 
form Supernovae that are intrinsically "rare" events. 
(iv) In the analysis of surveys, in special for the case 
of galaxies with emission lines produced by massive 
stars, that are intrinsically a small fraction of the 
total number of stars but dominate the emission.... 
In all these cases, the use of multiwavelength obser- 
vations and self-consistent synthesis models that in- 
cludes such statistical effects and computes the cor- 
relation between different wavelengths will be an as- 
set to constrain safely the physical properties of stel- 
lar systems (form small star forming regions to high- 
redshift galaxies). The current model outputs can 
be found in |rttp : //www . laef f . esa.es/mcs/SED/. 
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